Relationship between monomer packing and receptor binding domain pocket status in the spike trimer of SARS-CoV-2 variants (preprint)

Existence of a SARS-CoV-2 spike protein trimer form with closer packing between monomers when receptor binding domains (RBDs) are all down, locked as opposed to closed, has been associated with linoleic acid (LA) binding at neutral pH, or can occur at acidic pH in the absence of LA binding. The relationship between degree of closure of the LA binding pocket of the RBD, and monomer burial in the trimer, is examined for a range of spike protein structures, including those with D614G mutation, and that of the Delta variant (which also carries D614G). Some spike protein structures with this aspartic acid mutation show monomer packing approaching that of the locked form (at neutral pH, without LA binding) for two segments, a third (around the RBD) remains less closely packed. Analysis of other coronavirus RBD structures suggests that mutation of the RBD in spike protein of the Omicron variant could lead to LA binding pocket changes. It is proposed that these changes could lead to one of two consequences for the Omicron variant spike protein (which also has the D614G mutation), at neutral pH and without LA binding, either easier access to a locked form throughout that leads to cooperative transitions between all RBD down and all RBD up, or maintenance of a spike trimer with locked characteristics C-terminal to the RBD at the same time as the RBD is free to transit between down and up states. The situation may also be impacted by spike protein charge mutations in the Omicron lineage that alter pH-dependence around the RBD, in a similar way to the changes induced elsewhere by D614G.

Predicted pH-dependent stability of SARS-CoV-2 spike protein trimer from interfacial acidic groups (preprint)

Transition between receptor binding domain (RBD) up and down forms of the SARS-CoV-2 spike protein trimer is coupled to receptor binding and is one route by which variants can alter viral properties. It is becoming apparent that key roles in the transition are played by pH and a more compact closed form, termed locked. Calculations of pH-dependence are made for a large set of spike trimers, including locked form trimer structures that have recently become available. Several acidic sidechains become sufficiently buried in the locked form to give a predicted pH-dependence in the mild acidic range, with stabilisation of the locked form as pH reduces from 7.5 to 5, consistent with emerging characterisation by cryo-electron microscopy. The calculated pH effects in pre-fusion spike trimers are modulated mainly by aspartic acid residues, rather than the more familiar histidine role at mild acidic pH. These acidic sidechains are generally surface located and weakly interacting when not in a locked conformation. In this model, their replacement (perhaps with asparagine) would remove the pH-dependent destabilisation of locked spike trimer conformations, and increase their recovery at neutral pH. This would provide an alternative or supplement to the insertion of disulphide linkages for stabilising spike protein trimers, with potential relevance for vaccine design.

Uncovering cryptic pockets in the SARS-CoV-2 spike glycoprotein (preprint)

The recent global COVID-19 pandemic has prompted a rapid response in terms of vaccine and drug development targeting the viral pathogen, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In this work, we modelled a complete membrane-embedded SARS-CoV-2 spike (S) protein, the primary target of vaccine and therapeutics development, based on available structural data and known glycan content. We then used molecular dynamics (MD) simulations to study the system in the presence of benzene probes designed to enhance discovery of cryptic, potentially druggable pockets on the S protein surface. We uncovered a novel cryptic pocket with promising druggable properties located underneath the 617-628 loop, which was shown to be involved in the formation of S protein multimers on the viral surface. A marked multi-conformational behaviour of this loop in simulations was validated using hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments, supportive of opening and closing dynamics. Interestingly, the pocket is also the site of the D614G mutation, known to be important for SARS-CoV-2 fitness, and within close proximity to mutations in the novel SARS-CoV-2 strains B.1.1.7 and B.1.1.28, both of which are associated with increased transmissibility and severity of infection. The pocket was present in systems emulating both immature and mature glycosylation states, suggesting its druggability may not be dependent upon the stage of virus maturation. Overall, the predominantly hydrophobic nature of the cryptic pocket, its well conserved surface, and proximity to regions of functional relevance in viral assembly and fitness are all promising indicators of its potential for therapeutic targeting. Our method also successfully recapitulated hydrophobic pockets in the receptor binding domain and N-terminal domain associated with detergent or lipid binding in prior cryo-electron microscopy (cryo-EM) studies. Collectively, this work highlights the utility of the benzene mapping approach in uncovering potential druggable sites on the surface of SARS-CoV-2 targets.

A model for pH coupling of the SARS-CoV-2 spike protein open/closed equilibrium. (preprint)

SARS-CoV-2, causative agent of the COVID-19 pandemic, is thought to release its RNA genome at either the cell surface or within endosomes, the balance being dependent on spike protein stability, and the complement of receptors, co-receptors and proteases. To investigate possible mediators of pH-dependence, pKa calculations have been made on a set of structures for spike protein ectodomain and fragments from SARS-CoV-2 and other coronaviruses. Dominating a heat map of the aggregated predictions, 3 histidine residues in S2 are consistently predicted as destabilising in pre-fusion (all 3) and post-fusion (2 of 3) structures. Other predicted features include the more moderate energetics of surface salt-bridge interactions, and sidechain-mainchain interactions. Two aspartic acid residues in partially buried salt-bridges (D290 - R273 and R355 - D398) have pKas that are calculated to be elevated and destabilising in more open forms of the spike trimer. These aspartic acids are most stabilised in a tightly closed conformation that has been observed when linoleic acid is bound, and which also affects the interactions of D614. The D614G mutation is known to modulate the balance of closed to open trimer. It is suggested that D398 in particular contributes to a pH-dependence of the open/closed equilibrium, potentially coupled to the effects of linoleic acid binding and D614G mutation, and possibly also A570D mutation. These observations are discussed in the context of SARS-CoV-2 infection, mutagenesis studies, and other human coronaviruses.

protein-sol pKa: prediction of electrostatic frustration, with application to coronaviruses (preprint)

Evolution couples differences in ambient pH to biological function through protonatable groups, in particular those that switch from buried to exposed and alter protonation state in doing so. We present a tool focusing on structure-based discovery and display of these groups. Since prediction of buried group pKas is computationally intensive, solvent accessibility of ionisable groups is displayed, from which the user can iteratively select pKa calculation centers. Results are color-coded, with emphasis on buried groups. Utility is demonstrated with coronaviruses, which exhibit variable dependence on the acidic pH of the endocytotic pathway. After benchmarking with variants of murine hepatitis virus, a pair of conserved histidine residues are identified that are predicted to be electrostatically frustrated at acidic pH in a common structural core of pre- and post-fusion coronavirus spike proteins. We suggest that an intermediate expanded conformation at endosomal pH could relax the frustration, allowing histidine protonation, and facilitating conformational conversion. This tool is available at http://www.protein-sol.manchester.ac.uk/pka/.